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MATHEMATICAL  PROCESSING  OF  RANGE  AND  RANGE 
RATE  TRACKING  DATA 


Wang  Zhengming  Zhou  Haiyin 

ABSTRACT  Employing  radar  or  optical  equipment  to  carry  out 
tracking  measuremLts  with  regard  to  spacecraft,  it  is  possible  to 
obtain  observation  data  associated  with  ranges  and  their  rates  of 
chanae  Due  to  causes  in  many  areas  such  as  measurement  equipment 
oreSision,  environmental  conditions,  and  so  on,  these  measurement 
data  ii?lide  errors  in  all  cases.  Besides  random  errors  associated 
with  ranae  finding  and  range  change  rates,  there  are  also  systemic 
errors  associted  with  range  finding.  Through  the  setting  up  of 
appropriate  mathematical  models,  spline 

deany  siow  that 

the  methods  in  question  possess  very  high  accuracy. 

SUBJECT  TERMS  Spacecraft  tracking  Linear  system  Spline 
function  Error  analysis 


1  INTRODUCTION 


No  matter  what  type  of  measurement  system  is  used,  carrying 
out  tracking  measurements  with  regard  to  spacecraft  requires,  in 
all  cases,  measurements  [1]  of  ranges  and  their  rates  of  change. 
Because  of  many  types  of  complicated  reasons,  as  far  as  measurement 
values  are  concerned— besides  containing  random  errors— range 
finding  also  includes  system  errors [1-3].  These  errors  severely 
influence— even  to  the  point  of  distorting — the  true  image  of  the 
facts.  Also,  because  measurements  are  often  indirect  (for  example, 
MISTRAM  systems),  when  calculating  parameter  data  which  is 
necessary  for  use— during  transmission  processes  determined  by^ 
calculation  formulae — errors  are  often  enlarged  from  several  times 
to  several  tens  of  times.  This  is  not  permissible  in  a  number  of 
situations  where  there  is  a  need  for  high  precision  measurement 
data,  for  example,  those  in  such  things  as  missile  accuracy 
analyses,  spacecraft  orbital  predictions  and  so  on  [4]. 

With  regard  to  making  use  of  measurement  data  to  give  high 
precis ion  estimates  of  ranges  and  their  rates  of  change,  systemic 
errors  associated  with  estimated  range  finding  are  problems  jointly 
concerned  with  measurement  units  (precisions  of  equipment 
concerned)  and  data  utilization  units  (concerned  with  accuracies  of 
range  measurements  and  their  rates  of  change). 

This  article  considers  problems  in  launch  coordinate  systems • 
yiQ  first  of  all,  demonstrate  that — with  respect  to  spacecraft  in 
stable  flight  configurations  (no  stage  separations)  observed 
ranges  and  their  rates  of  change  are  capable,  in  all  cases,  of 
being  expressed  by  the  use  of  a  cubic  spline  and  its  derived 
functions.  On  the  basis  of  system  error  models — making  use  of 
spline  function  theory  and  regression  analysis  methods  and  through 
the  setting  up  an  appropriate  mathematical  model  mathematical 
methods  are  given  which  are  not  only  able  to  estimate  ranges  and 
their  rates  of  change  but  are  also  able  to  ^  estimate  range  finding 
system  errors.  Theoretical  analysis  and  simulation  calculations 
both  clearly  show  that  the  methods  in  question  have  trimmed  errors 
which  are  small,  and  calculations  are  very  convenient. 

2  SPLINE  FUNCTION  DESCRIPTIONS  ASSOCIATED  WITH  RANGES  AND  THEIR 
RATES  OF  CHANGE 


Assuming  oxyz  to  be  the  launch  coordinate  system,  spacecraft 
orbit  parameters  are  x(t),  y(t),  z(t), 

Here,  x(t),  y(t),  z(t)  express  the  position  of  spacecraft  at 


instant  t. 


is  the  speed  of  spacecraft  at 


instant  t  .  This  article  adds  dots  to  the  top  of  dependent  /18 
variables  to  stand  for  the  derivatives  of  the  functions  in  question 
with  respect  to  time.  The  number  of  dots  on  top  stands  for  the 
order  of  the  derivative. 


X 


Assuming  that  the  coordinates  of  observation  stations  in 
launch  coordinate  systems  are  (xo/yo/zo)^  then/  observed  ranges  ares 

jj(t)  =  )*+(»-»,)*+ (z-z,)*  ( 1 ) 


Assuming  that  xt/y(t)/Z(t)  are  all  fourth  order  continuous 
differentiable  functions,  then,  from  equation  (1),  it  is  possible 
to  deduce  that  R(t)  is  also  a  fourth  order  differentiable  function. 
Making  use  of  x(t) /y(t) /Z(t)  as  well  as  the  first,  second,  and 
third  order  derivatives  of  special  points,  direct  solution  is  made 
with  respect  to  R(t).  It  is  possible  to  demonstrate  that  R(4)(t) 
is  a  very  small  number. 

Below,  we  study  data  processing  problems  associated  with 
sections  where  time  length  is  [To,Tm].  Using  subscripts  To=tl, 
Tm=tn,  T  j=To+h*  j ,  j=-l ,  0 . . .  ,m+l .  ||  ||  •  =max  . 


Using  Sj(t)  to  represent  the  jth  spline,  Bj(t)  records  the  jth 

standard  B  spline.  . 

Lemma  1  [5]  Assuming  R(t)€C4[To,Tm]  and  S(t)  satisfies  the 

conditions 

I J?(0  -  S,  it)  I  II  it)  II  --AS  m)-  .9,(01  <  ^  II  R^*\t)  II  ( 2 ) 

of  cubic  interpolative  spline  functions,  one  then  has 


The  first  order  and  third  order  standard  B  splines  which  this 
article  uses  are  respectively  defined  as: 


W  >1 

|t|,  M<1- 


B,(t) 


0,  W  >2, 

m<i 

2  3 

— i^+T*-2M 


1<  M  <2 


Lemma  2  There  exists  a  unique  set  of  coefficients  (b-1,  bo, 
bl , . . . ,  bm+1 )  so  that 
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(3) 


n+  I 


S,(t)=  E  b,By 


h 


satisfies^condi^^  the  two  lemma  above,  it  is  possible  to  obtain 
Theorem  1  During  the  normal  flight  phase  of  spacecraft 
(refers  to  phase  with  no  stage  separation),  tracking  range  R(t)  of 

spacecraft  and  the  rate  of  change  R(0  can  be  described  by  the 

use  of  cubic  spline  function  S3(t)  as  well  as  the  derivative  S,(0 
So  Iona  as  m  is  sufficiently  large,  making  h  (tn-tl)/m 
sufficiently  small,  it  is  then  possible  to  guarantee  trimmed 
being  sufficiently  small.  In  this,  S3(t)  is  given  by  equation  (3). 


3  MATHEMATICAL  MODELS  OF  OBSERVATION  DATA 

Observation  data  associated  with  R(t)  can  be  represented  as 


(4) 


In  this  R(ti)  is  the  true  value  of  R(t)  for  instant  ti  ;  C(ti)  is 
JSe  observation  system  error  associated  with  range  finding  at 
instant  ti  ;  t^  is  random  error. 

Mt)  observation  data  can  be  represented  as 

p,  =  Mti}+5,  (5) 

In 

this  MU)  is  the  true  value  of  at  instant  ti; 

Ti  ls  rando^  error  associated  with  ® C(t) 

Below,  we  take  problems  associated  with  R(t),  Mt)  /  and  C(t) 
and  turn  them  into  parameter  estimation  problems. 

As  far  as  the  consideration  of  problems  within  the  time  phase 
t?o  Is  cSncernsd,  assuming  that  ti  is  the  observation  instant, 

k  is  the  number  of  observation  iterations  each  second,  and 

1- i=To+^ i”l ) /k,  i=l,2,...,n,  then,  ji 

ti  TO  U  }/  ,  Tj=tkhj+1,  j=0,l,...,m;  n=khm+l 

On  the  basis  of  engineering  background  (see  references  [1- 

3)),  within  [To,Tm],  C(t)  can  be  described  by  the  use  of  first 

order  spline  number  Sl(t)i 


(6) 
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From  discussions  in  2,  we  already  have 


■  ^1 


=  b,Bl 

»--i  ' 


t-Ts 


)■ 


ii(0> 


M4>1 

!i,(o  =  E 


B. 


(^) 


(7) 


Because  of  this,  from  equations  (4)  ~  (7), 

to  obtain 


it  is  possible 


(8) 


If  we  are  able  to  obtain  from  equation  (8) 
estimates  for  (b-l,bo, . . . ,bm+l )  as  well  as  (do,dl, . . . ,dm) ,  then, 
we 

are  able  to  respectively  obtain  estimates  for  C(t),  R(t),  and  R(t) 
from  equation  ( 6 )  and  equation  ( 7 ) . 


4  ESTIMATION  METHODS  AND  ERROR  ANALYSIS 

It  is  clearly  shown  by  theoretical  analyses  and 
simulation  calculations  that,  during  relatively  stable  flight 
phases,  h  can  be  selected  as  5  ~  10s.  It  is  possible  to  guarantee 
that  II  5,«)  II  .  ,11  II  -  ,11  S,«)-C(0  |I  ., 

are  all  very  small  (primarily  comparatively  random  errors  as  well 
as  being  within  the  ranges  permited  by  engineering  background) . 
During  orbital  tracking,  each  second,  it  is  possible  to  obtain  over 
20  observed  data.  Thus,  n^lOOm. 

On  the  basis  of  the  special  characteristics  of  measurement 
equipment,  in  relation  to  random  errors,  it  is  possible  to  assume 
that 


Fe,  =0,F5,  — 0,£'f(5,  —  0,  I,/  — l,2,*"',n  i  (9) 

£ef=a|,Fa|=0|,EC|C,=O,E5,a,=O  f#;  J 

Here,  ai,  0i  (i=l,2, . . . ,n)  are  already  known  positive  numbers.  It 
is  possible  to  obtain  them  from  appraisals  of  measurement  equipment 
precisions . 


Below,  use  is  made  of  observation  data  and  model  (8)  ot  give 
estimates  of  parameters  (do,  dl,...,  dm)  as  well  as  (b-1,  bo,..., 
bm+1) . 

Lemma  3  [5]  Assume  A  is  any  constant.  Then,  A  can  be 
expressed  as: 

^='p_AB,{l=p-yt  .ecT..r.3 


From  Lemma  3,  it  is  possible  to  know  that  parameters  which 
await  estimation  in  model  (8)  equations  are  not  completely 
independent.  We  take  model  (8)  and  rewrite  it  as  /20 


(10) 


Here,  use  is  made  of  Lemma  2  as  well  as 


Note 


a,:=b,  +  dt,  — ,>n+l 

gj  =  dj-A,  ;  =  l,2,-,m 

X  —  (j«>  » (i«+j) 


,m  +  3 


JLb,  «-=U2,-,«,  ;=1.2,-, 

01  \  n  /  , 

,  _  I  1  Ip  ~ *V  i  =  n+l,?i+2,*'’»2n,  ;  = 

fcT  *V  h  r 

1  p  ,'=;i,2,  — ,n,  ;=m  +  4,-,2m+3 

at  \  k  ' 

\0,  i  =  «  +  l,  — ,2n,  ;=m  +  4,  — ,2m  +  3 

5  =  *"  f  ^■+1 

v_r».  y-  ^V, 


771+3 
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in  this  way,  we  obtain  linear  regression  model: 

r=X/?  +  c  I 

£c=:0,COY(e)-7,,Ki.  ^ 


(11) 


Theorem  2  Assume  that  random  errors 
finding  and  its  rate  of  change  cl  and  61 
conditions  of  equation  (9).  Then,  it  is 
regression  model  (11).  Moreover,  square 
defined  matix.  Thus,  from  formula 

^ts  =  iX'Xy^X'Y 


associated  with  range 
respectively  satisfy  the 
possible  to  obtain  linear 
matrix  X’X  is  a  positively 

(12) 


it  is  possible  to  give  linear  unbiased  estimates  associated  with  {} 

variances  which  are  uniformly  minimal.  _  vnnw  i-hai-  the 

Proof  From  equation  (9),  it  is  possible  to  know  that  the 

second  form  of  equation  (11)  is  set  up.  From  discussions  in  2, 
?Sed  e??ors  associated  with  using  S3(t)  to  appr ox^ate  R ( t ) , 
usina  SiU)  to  approximate  li(0  and  using  Sl(t)  to 
approximate  C(t)  can  be  ignored  in  calculations.  Thus,  the  fi 

form  in  elation  then,  by  the  Gaues- 

estimates  consitent  with  parameter  i?  variances  associated  wit 

linear  regression  model  (11).  vx  is  oostivelv 

Thus,  we  only  need  to  prove  that  matrix  X  X  is  posriveiy 

defined.  Note 

\_(X,  Xa 


In  this. 


V  Tft  a  ^ 2Tn+4 ^  row  by  (2ni+3)  column  matrix  formed 

f^inhl'first^row  ofx?  ^hftt/2  rowT  ihe  hij  row  (j=l,2 . m), 

the  mkh+1  row,  the  mkh+(kh/2)  row,  and  the  kn(m+j)  ,  - 

j;-ta-io^i;e.  a/'the 

"ea«e<?riorairAt  tf’ 

IZX  llZl  is:?:  e?ln;gard  to  any  number  net  r  . 
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/ 


X  r . 


all  have  /X’Xr  >( illegible )  /  X  ' 


far  as  another  definition  from  X  is  concerned/  we  have; 


(  X 

6 

2 

T 

1 

T 

0 

0  ••• 

0 

0 

0 

1 

48 

23 

48 

23 

48 

I 

48 

0  .  - 

0 

0 

0 

0 

1 

6 

2 

3 

1 

6 

0  ••• 

0 

0 

0 

0 

0 

1 

6 

2 

T 

i  ... 

6 

•  • 

0 

0 

0 

• 

• 

? 

• 

T*  .  * 

i 

2 

1 

0 

0 

.0 

0  ••• 

6 

T 

by 

/O  0 

i « 
1  0 
0  1 


^0  0 


0^ 

0 

0 

0 


-JL  0  \  0 

Zk  2h 

_  j _ L  -L  JL 

16A  8A  8A  16 


0 

0 


— L  0  i  0 

Zh  2k 


0  0 
0  0 
0  0 


I  0 


0  0  0 


0  J- 

Zk  2h) 


Making  use  of  elementary  transformation  methods,  it  is  possible  to 
p?oie\hat  X  is  a  full  column  ordered  matrix.  Thus,  X 
X  is  positively  defined.  Because  X  X  is  positively 

defined,  the  theorems  are  proven. 

Below,  we  come  to  discussion  of  estimation  errors  associated 
with  R(t)+C(To),  A(t)  /  as  well  as  C(T)-C(To).  Assuming 

J?(0+C(r,)  =  S,«)  =  I] 

I  -  - 1 

C(o-C(r.)  =  s,«)  =  Eg>B,(i=^) 

/.I  •  ' 

note  /9m  =  (  — ^ 


^  /  l—T,  \ 

Use  S,(0  =  E  aiH-r-  h 

i  —  1  \  n  / 

to  act,  respectively,  as  estimates  of  S3(t), 


Si(t) 


,  and  Sl(t). 


Define 


£(S,)=£i]  E(Sty=EYl  |i's(<i)-i’s(<«)! *, 

f(S,)  =  £Z]  |S,(t,)-S,(^)I  *,  a  =  g-(gi,gzr"tSny  1 


Z=(Zi.  i),x  (■+j).  ^  =  />•«<•+».  tU-(Ui,  ,).«■. 


In  this, 


Because  of  this,  ,1 j(s.)iE||  R’c'S-a)  ||  ■,  E(S.)  =E  II  £/(,-,)  T 

Leituna  4  [6]  Under  the  assumptions  of  model  (11),  one  has 

xr\  cov(a)=r.,  cov(J‘)=7/, E(/s«)_r=/?,  EGy=a,EG'>=g>  COY  (pLs)  =a' 

H^hre  rx’XWl  is  a  2m+3  order  square  matrix.  VI  is  an  m+3  order 
matrix.  It  is  composed  of  elements  at  locations  where  (X  X)-l 
front  m+3  lines  and  forward  m+3  columns  intersect.  V2  is  an  m 
order  matrix.  It  is  composed  of  elements  at  locations  where 
(X'X)-l  rear  m  lines  and  rear  m  columns  instersect. 

Theorem  3  Under  the  assumptions  of  model  (11),  we  have 
S3)=tr(Z  '  ZVl) ,  E(  Sj  )=tr(W'WVl),  and  E(Sl)=tr(U'UV) .  Here,  .<r(*> 
stands  for  the  sum  of  the  elements  of  primary  diagonals  of  square 

Proof  Only  proof  is  of  the  expression  for  E(S3).  The  rest 
are  similar. 


E(St)=E  \\  Z(.7-a)  II  *  =  Eia-ayZ'Z  {a-a) 


-EtKU-a)'Z'Zia-ay^  =  EirCZ'Z  (fl-a)'( 

=<r  fZ'ZCOV(  «))=<r(Z'Z7,) 

What  is  worth  pointing  out  is  that  the  conclusions  described 
above  are  obtained  under  the  presupposition  that  trimmed  errors  can 
be  assumed  to  be  capable  of  being  ignored  in  calculations.  As  far 
as  how  to  guarantee  this  point  is  concerned,  it  is  possible  to  make 
use  of  simulation  methods — simulating  the  production  of  an  orbit 
(for  example,  using  actual  orbit  design  tracks)  as  well  as  system 


9 


errors.  Assuming  C(t)=at(a  is  a  constant),  calculations  are  done 
of  R(t).  For  cases  where  random  errors  are  not  added,  calculations 
are  done  of  S3(t)  and  Sl(t).  In  conjunction  with  this,  solutions 
are  gotten  for 


i-i  .-1  appropriately 

selecting  h  to  guarantee  the  above  three  summations  are, 
respectively,  smaller  than 

ifo  S 


Through  large  amounts  of  simulation  calculations — during  stable 
flight  phases — selecting  h=5s,  it  is  possible  to  guarantee  the 
requirements  described  above  being  satisfied. 


5  SIMULATION  CALCULATION  EXAMPLES 


Example  1 


Assuming  that  orbital 
x(i),  yit),  z(r) 

equations 


pareuneters  *(<),  Vit),  ar(0,x(t), 
satisfy  orbital  movement 


Au 

At 


+8f  +«., 


(13) 


In  these,  n  is  one  of  x(t),  y(t),  or  z(t). 

are,  respectively,  resistance  acceleration,  acceleration  of 
gravity,  drag  acceleration,  and  Ke  acceleration. 

Given  initial  values  associated  with  the  instant  To  [x(To,..., 
2.  (To)]',  from  equation  (13),  it  is  possible  to  produce  /23 

orbits.  Thus,  it  is  possible  to  calculate  true  yalues  R(ti)  and 
D  (ti)  associated  with  measured  data  for  various  instants  ti. 
Assuming  that  range  system  error  is  C(t)=at — in  situations 
where  random  errors  are  not  added — one  obtains  observation  data  yi 
=  R(ti)  +  Mti,  pi  =  R(t,),.  Table  1  below  gives  the  differences 

between  estimated  values  and  real  values  obtained  by  the  use  of  the 
methods  in  this  article. 

Trimmed  error  square  summations  are,  respectively. 


Table  1 


Trimmed  Errors  Displayed  by  R(t),  C(t),  and  R(t)Splines 


i 

R(ti)— R(ti) 

1 

7.574  85D- 5 

101 

5.172  02D-3 

201 

6.442  70D-3 

301 

4.707  85D-3 

401 

•  5.650  41D- 3 

501 

4,512  03D-3 

601 

—1.033  14D-4 

c  (/i)-cai) 

0 

-2.570  15D-3 
-3.967  80D-3 
-2.830  ISO- 3 
-5.002  75D-3 
-4.405  23D-3 
-2.304  72D-3 


5.305  57D-3 
-3.166  62D-4 
1.282  69D-4 
-3.003  38D-5 
3.677  13D-5 
3.422  080 -4 
-4.240  71D-3 


Examole  2  On  the  foundation  of  Example  1, 

produSHcti)  and  f  (ti)  with  random  errors  added  to  obtain 

observation  values,  yi  -  R(ti)  +C(^)  +  €i,  pt—itUi)+dt, 

^  Kin  o  0121*  5,~iV(0.0.03*).032);  i  Ee,ej-Q,  E6i6j  =  O, 

In  these,  ni  As  related  to  random  error 

ai^mptioAr^it  is  piLibli'io'coAsult  reference  [1]. 

Table  2. 

Table  2  Estimation  Errors  Associated  with  the  Methods  of  This 
Article 


We  obtained 


1 

101 

201 

301 

401 

501 

601 


1.877  43D- 2 
2'155  90D-2 
3.293  63D-2 
3.059  65D- 2 
2.539  76D-2 
8.181 96D- 3 
-4,683  51D-3 


-2.040  03D-2 
-5.126  47D-2 
- 1.426  44D-?. 
-  1.948  12D-2 
-i.262  67D-3 
-2.238  30D-2 


-8.308  79D-3 
-3.992  08D-3 
5.864  7  D  -4 
-4.633  22D-4 
-2.134  26D-3 
-3.871  58D-3 
1.89023D-4 


Estimation  error  square  summations  are,  respectively, 

*  =  0.284  950,  X]  = 

i  - 1 


2  *  =  1.031  883D-2 


i*  1 


The  methods  of  this  article  are  not  only  suitable  for  use  in 
forms  where  system 

assumed  that  ranp  rs  then  the  methods  of  this 

artil!fare  approp^ia^ffor  use  in  “y  °th4r  situation^associated 

triaL"s;i"eiSnclionil(trannhange  it  into  an  /24 


II 


appropriately  high  order  spline  function  in  order  to  guarantee 
I II  ^(1)  —  5j(t)  II being  sufficiently  small  and  est^ation  precisions 
not  being  influenced.  Realizations  of  theoretical  analyses  and 
actual  calculations  are  similar. 
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